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ABSTRACT 


Fully resolved intensity profiles of various lines in the C0 9 bands 
at 9. A p m and 10. A M m have been measured on Mars using an infrared 
heterodyne spectrometer. Analysis of the line shapes shows that the Mars 
atmosphere exhibits positive gain on these lines, providing the first 
definite detection of natural optical gain amplification and enabling 
identification of these lines as the first definite natural laser ever 
discovered . 
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Since their Invention some 25 years ago, optical lasers have become 
nearly ubiquitous tools in the laboratory and in everyday life (e.g. 
super-market price scanners). Despite this somewhat deceiving abundance 
of man-made lasers, no definite naturally-occurring lasers have ever been 
reported, even though the extreme variety of physical and chemical 
environments of extra-terrestrial objects plausibly argues that natural 
lasers must exist (1). Natural microwave amplifiers (masers) have been 
found to be abundant in interstellar clouds and some circumstellar 
shells, primarily among the rotational level populations of certain 
molecules (e.g. OH, SiO, H^O), and are all characterized by the property 
hv < kT, where T is the kinetic temperature. However, optical lasers 
are characteristic of electronic or vibrational transitions for which 
hv >> kT. Many examples of natural non-thermal optical emission have 
been found, such as the infrared and ultraviolet aurorae and/or day glows 
of Earth, Jupiter, Mars, and Venus, but it has never been established 
that a population inversion exists for any of these emissions. The 
relative populations of the two molecular levels must be inverted in 
order for gain amplification to occur, the essence of lasing. V/e report 
here the discovery of a population inversion and natural gain 
amplification in the 10. 4 um CO^ laser bands on Mars, representing to our 
knowledge the first definite identification of a natural infrared laser. 

The observations were made with the Goddard infrared heterodyne 
spectrometer (2) during the period January-April , 1980 when the planet 
was near opposition. The beam size (HPBW) was 1.7 arc-sec and the disk 
of Mars was «/> 13.8 are-sec diameter, providing good spatial resolution on 
the planet. Fully resolved atmospheric line profiles were measured at 
various locations on the disk, including many individual 
rotational-vibrational lines in both the 10.4 um (00°1 - [10°0, 02°0]^.) 
and 9.^ u m (00°1 - [10°0, 02°0] TT ) bands. Analysis of the full set of 
lines is proceeding, and we report here some results obtained from 
analysis of several measurements of the 10.33 u™ R8 line at 967.7072cm - " 1 . 

The intensity profiles were measured simultaneously at 25 MHz 
(0.0008 cm -1 ) resolution with 64 consecutive channels, and at 5 MHz 
(0.00016 cm -1 ) resolution with a second bank of 64 consecutive channels. 
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It is important to note that all 128 channels were recorded 
simultaneously, thereby eliminating registration errors and drift. 
Absolute Intensity calibration and removal of terrestrial atmospheric 
lines was achieved by lunar comparison spectra, measured nearly 
simultaneously and scaled to the appropriate air mass. 

Fully resolved intensity profiles of the R8 line are show) (Fig. 1) 
for beams placed at disk-center, the north and south polar regions, and 
the equatorial east and west limbs (3). ileveral qualitative conclusions 
are evident by inspection of toe spectra. First, In toe absence of 
significant continuum opacity (local dust storms or clouds), the emergent 
intensity far from line center is a direct measure of thermal emission 
from the surface and therefore of the surface temperature. The measured 
surface temperatures agree well with data obtained by the Mariner 9 and 
Viking orbiters, where comparisons have been made to date. Second, 
analysis of the broad wings of the absorption line enables extraction of 
various atmospheric parameters by inversion of the radiative transfer 
equation for each location. Third, the bright emission at line center is 
evident at all positions with an intensity which is variable with 
location (U). It is this core emission which we will show constitutes a 
natural infrared laser. 

Quantitative analysis of the disk-center spectra was performed in 
the following way. We first analyzed the line wings, from the continuum 
to within 50 MHz of line center. The solid surface temperature 
(typically 250K) was determined from the absolutely calibrated 
intensities far from line center. A discontinuity of 30K between the 
solid surface temperature and the gas temperature near the surface was 
assumed, in accord with radiative equilibrium and spacecraft 
observations (5). A mid-latitude model temperature profile typical of 
local noon, as determined by spacecraft measurements, was adopted and is 
shown in Fig. ?. A temperature lapse rate of 2.0 K/km was assumed to an 
altitude of ^5 km, followed by an isothermal atmosphere (typically 150K) 
above «5 km. The shaded portion of Fig. ?. represents the range of 
temperatures measured by the Viking probes in the so-called isothermal 
region (6) and the mean lapse rate is also in agreement with measurements 
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by Mars 6 , and Viking for local noon (5). 

We initially assumed that the (00°1) state was in 
local-thermodynamlc-equlllbrium throughout the lower atmosphere (2 < 

45 km) and performed an Iterative solution to the radiative transfer 
equation as a function of surface pressure. Two lines were well-fitted 
with these assumptions (e.g. Fig. 3a), however the remaining three lines 
showed too little atmospheric self-emission within *300 MHz of line 
center. Mo reasonable set of lapse rates, surface temperatures and 
pressures, or haze opacity produced a suitable fit for these three lines. 
However, all three remaining lines are well-fitted by the above described 
model if the (00°1) source r f iction is allowed to fall below the LTE 
value in the lower atmosphere of Mars. Furthermore, the vibrational 

_ 1 C -? 1 

relaxation rate coefficient derived in this way (K * *4.7x10 em'sec ) 
agrees well with laboratory values (7). In addition, the retrieved 
surface pressures for all five lines agree well with the known values 
determined by the Mariner 9 topographic maps and the Viking Lander annual 
pressure measurements. We therefore believe that our observations 
represent the first observation of the failure of LTE in a planetary 
tropospnere and we attribute this to screening of radiative exchange in 
the urn band by airborne dust. Radiative energy exchange in a clear 
CO^ atmosphere would keep the (00°1) level in LTE to an altitude many 
scale heights above the level at which the collisional and radiative 
lifetimes were comparable ( 8 ), however a small amount of dust will 
efficiently capture quanta from the optically thick *4.3 urn band, thereby 
screening this exchange. T t is beyond the scope of this paper to discuss 
this effect further, except to note that the detailed physics occurring 
in the troposphere does not significantly affect the principal conclusion 
of this paper. It is sufficient to note that we have successfully 
modelled the emergent tropospheric intensity, enabling us to strip this 
from the total emergent intensity, leaving the residual core emission. 

Two fitted examples of the main lines are shown in Fig. 3 . Tn 
addition to the calibrated single side-band experimental data we show 
separately the modelled self-emission of the lower atmosphere, the 
modelled transmitted surface intensity, and the sum i.e. the modelled 
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emergent intensity. We next stripped the modelled emergent Intensity 
from the observed and fitted the residual core emission with • Gaussian 
profile, using a least-squares analysis for the amplitude, width, and 
velocity shift. Typical results are shown In Fig. 9, 

The derived core emission line-widths are a direct measure of the 

kinetic temperature in the emitting region, and the integrated line 

brightness (E.) can be related directly to the column density In the 
” _ 

upper state (00 1) by the following relation: 
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where J' = 9, the rotational quantum number of the upper level of the R8 
line. We assune the core emission to be emitted isotropically, which 
will be justified a-posteriori if the deduced atmospheric gain is 
sufficiently small that directional effects within the emitting region 
may be neglected. The effect of this assumption will be to underestimate 
the total atmospheric gain. We take the kinetic and rotational 
temperatures to be equal in the emitting region which has been shown to 
be peaked at «r75 km (9). The radiative transition rates (A^, A^ 0 > for 
the 9.9 nm and 10.9 u m bands have been measured in the laboratory (9). 

The rotational partition function may be taken to be 0.69 T/B to 

good approximation. The results of this analysis are given in Table I 
for five observations of the RS line at disk center. We find that the 

column density in the (00°1) state is not less than 2.0X10 1 ** cm~^ ± 

13 -2 

2.0x10 cm '. 

Uc note that one line (Chryse Planitin) shows a substantially higher 
kinetic temperature (15 IK) than the remaining four observations which 
have a mean kinetic temperature of 1 1 3 ± 10K . It’s width (19.25 MHz) is 
nearly three standard deviations greater than the mean width of the other 
four lines ( 16.5+0.8). We have carried through parallel analyses based 
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on four and five lines respectively .because the Viking entry probes 

measured considerable thermal structure in the region 60-120 km 

suggesting that the observed high temperature (15 IK) may be real. We 

shall show that inclusion of this line in the analysis does not affect 

the main conclusion of this paper. We note that the total energy 

- 2-1 

radiated in the 9 and 10 pm bands is 18.5 ergs cm s which agrees well 

with earlier results (4). It is interesting to note that the energy 
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radiated in the R8 line is </• 10 times greater than would be radiated by 
this region of the atmosphere were it in local thermodynamic equilibrium. 


The total directional gain, f(v), ong the line of sight may be 
calculated (10) from the relation 


A . 

r(v) = — A (N-,- N t „ Jtt ) g(v) 

8n 10 0 * 


(2) 


where the w's are state statistical weights, and g(v) is the line shape 
function. The emergent intensity will then be given by (for unsaturated 
amplified spontaneous emission) by 


I(v = n(exp r(v) - 1) 


where n is some constant. If N > N w'/w", the state populations are 

J J 

said to be inverted and the gain constant is positive definite, a 
necessary and sufficient condition for gain amplification (i.e. lasing) 
to occur. 

We evaluate the maximum possible column density in the lower state 
[ 10°0, 02°0 in the following way. We must establish three physical 
parameters, vi 2 t (1) the base altitude level above which the emission 
originates, (2) the rotational temperature of the lower (I) level, and 
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(3) the vibrational temperature of the I level. The base altitude level 
may be established either by determining the region where the pumping 
energy is deposited, or more severely by requiring that the observed 
kinetic temperatures be greater than the condensation temperature. 
Johnson et al (4) have calculated the former, finding that </50g of the 
core emission originates at altitudes greater than 75 km. We have 
independently verified their result. Condensation temperatures of 130 K 
and 113 K are reached at altitudes of 48 Ion and 65 km, respectively (6). 
We will show that a population Inversion exists for each lower altitude 
bound, and therefore that the exact choice of altitude bound is 
unimportant for our main conclusion. 


The vibrational exchange rate for I with CO (00 0) is nearly gas 
kinetic (11) and the radiative relaxation rate is * 3 see . Therefore, 
collisional relaxation and radiative relaxation become comparable at ^120 
km, far above the emitting region, and we expect collisions to dominate 
radiative effects at lower altitudes. We therefore take Tyjg = ^ROT * 
i.e. the I state must be in local thermodynamic equilibrium, life 
have already directly measured the kinetic temperature in the emitting 
region (Table I), and we may now calculate the maximum lower state (T) 
column density above our lower altitude bound for an exponential 
atmosphere. We find that in every case the lower state column density is 
less than the upper state column density (Table II), i.e. a population 
inversion exists. 


V/e may now calculate the (minimum) directional gain at line center 

from: 
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where J" ■ 8, J* * 9, and the line-shape function at line center ia 



and we take B 00°i * P I * °*39 ew ” 1 02). 

The results are given in Table II. Note that when using the source 
function altitude bound, we have halved the upper state column density 
before calculating the gain. The gain is positive definite in all cases 
(13) and for the most reasonable lower bound (source function), the 
population inversion is greater tha>: 50:1. 

It is important to note that in every step we have routinely made 
assumptions t*iich would minimize the total gain, nevertheless the derived 
gain Is positive to an accuracy of ten standard deviations, based 
primarily on the accuracy with which the upper state population is 
determined (±101). The derived gain would result in gain narrowing of 
the line by *M0 KHz, an unobservable amount compared with the kinetic 
line widths. However, the atmospheric gain should be substantially 
greater at large zenith angles than in the zenith direction and a search 
for macroscopic gain narrowing seems warranted. 

In conclusion, natural gain amplification has been observed in the 
mesosphere of Mars, representing to our knowledge the first definite 
identification of a natural infrared laser. 
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1 Pertiel Global Map of the CO^ R8 Line on Mere. Small cirelea 
represent the plaoement and relative sice of the instrumental 
field of view. Spectra are displayed with a linear intensity 
soale and 25 MHz (0,0008 cm"" 1 ) spectral resolution. Note 
particularly the globally variable surface (continuum) 
temperatures, the (tropospheric) broad absorption line, and 
the (mesospheric) bright emission core. For clarity, we have 
omitted the higher resolution data near line center. 

2 The model temperature profile used to analyse the bread 
absorption wings contributed mainly by the lower atmosphere. 
The altitudes at which the self-emission source function peaks 
for given frequency difference from line center are given on 
the right. Also shown are the two principal regions where the 
(00°1) state Is found to depart from local thermodynamic 
equilibrium (LTE). Pj denotes' the level for which the 
normalized source function falls to 0.5, denotes the level 
for which the source function abruptly returns to LTE. The 
laser emission source function peaks at * 75 km. 

3a R8 line profile near local noon over Chryse Planitia (CML 

41°). The model requires LTE throughout the lower atmosphere 
to fit the emergent Intensities from the continuum to within 
50 MHz of line center, where the strong laser emission is 
seen. ( ,♦ are 25 MHz; X are 5 MHz channels) data used to 

model the wing are marked as . The intensity scale is 

linear, relative to a 260K single side-band black body. 

3b R8 line profile near local noon, north-west of Tharsls (CML 
120°). A non-LTE model *s required in the lower atmosphere to 
fit the line. An abrupt return to LTE at the 2.4 mbar level 
is found. In this case, the source function falls below 0.5 
at 3.05 mbars. Intensity scale is linear, relative to a 266K 
single side-band black body. 
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Figure Ua Laser ami salon at the core of tha R8 Una ovar Chryse 
Planltla. Tha underlying curve Is tha total em e r g ent 
Intensity in tha absence of laser emission, modelled as 
described in tha text. Tha difference between the ub served 
and modelled intensities is fitted to a gausslan as discussed 
in the text. Jtt.se shown are the half-width at half-sax inuai, 
the derived kinetic temperature, and the Integrated energy in 
the laser line. Intensity scale is linear, relative to a 
single aide-band 26 OK black body. 


Figure hb Laser anlsslon north-west of Thsrsis Montes. See text for 

discussion of analysis. Intensity scale is linear, relative 
to a single side-band 266K black body. 
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CORE EMISSION ANALYSIS 


BEAN CENTER 
LONGITUDE, LATITUDE 

Observed 

HALF-WIDTH 

IHALF-MAX 

(MHZ) 

APPARENT 

KINETIC 

TEMPERATURE 

(K) 

OBSERVED . 

LINE ENERGY*' 
ergs am" s sr 

MINIMIM 

REQUIRED COLUMN 
DENSITY IN CO. 
00 1 STATE * 
CM"* 

SE of Acidalia 
Planitia : 2+24° 

16.81 

115.3 

3-8K-2) 

2.19(14) 

Chryse 

Planitia:4l°+23° 

19.25 

151.1 

4.46(-2) 

3.03(14) 

Tempo 

Fossae: 82 °«-25° 

16.81 

115.3 

3.44(-2) 

1.97(14) 

NW of Thar sis 
Montes: 1 20+22° 

15.51 

98.1 

3.93(-2) 

2.07(14) 

Amazon is 

Planitia: 158°+27° 

17.39 

123.4 

2.87(-2) 

1.71(14) 

mean of 5: 

17. 1511.36 

b) 120.0H9 

3 .70( — 2 ) 
l5.9(-3) 

c) 2. 17(14) 
15.4(13) 

mean of 4: 

16. 5310. 79 

b) H3H0 

3.5K-2) 

i4.8(-3) 

1.99(14) 

12.0(13) 


a) Corrected to the sub-solar point. 

b) The condensation temperature (CO^) is *11 OK at 70 Km. 
Viking temperature extremes (60-116 Km) were 


VI 1 33.6K (92 Km) ♦ 154.6 K (64 Km) 4:13 PM LMT 

V2 115 K (116 Km) ♦ 157 K (88 Km) 9:49 AM LMT 

•2 — 1 

c) The integrated radiance in the 9 and 10pm bands is 18.5 ergs cm s at 
the sub- solar point. 


COLONS DENSITIES 
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t k 
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tNggQ.CB ] 

UPPER . INVERSION 

^001’°* ^ ^ 

CENTER 

Condensation 

48 Kb 

120K 

< 1.1 X 

lO 1 * 

> 2.2 x 10 1 * 

♦1.1(14) 

♦1.45(-3) 

ft 

65 Kb 

113K 

< 4.0 x 

10 12 

> 2.0 x lO 14 

♦2.0(14) 

♦ 2.7K-3) 

Source Function 

75 Kb 

120K 

< 3.2 x 

10 12 

14 

> 1.1 X 10 

+1.1(14) 

♦ 1.38(-3) 

n 

75 Km 

113K 

< 1.1 X 

10 12 

> 1.01 x lO 1 * 

+1.0(14) 

+1 .37 ( —3 ) 
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CONCLUSION: THE ATMOSPHERE EXHIBITS POSITIVE GAIN FOR ALL CASES. 
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FIGURE 3a 
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